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Abstract 



On 

O . We compute the distribution of quasideuterons in doubly closed shell nuclei. 

The ground states of 16 and 40 Ca are described in Is coupling using a re- 

m 

alistic hamiltonian including the Argonne v f 8 and the Urbana IX models of 
two- and three-nucleon potentials, respectively. The nuclear wave function 
contains central and tensor correlations, and correlated basis functions the- 



ory is used to evaluate the distribution of neutron-proton pairs, having the 
deuteron quantum numbers, as a function of their total momentum. By com- 
puting the number of deuteron-like pairs we are able to extract the Levinger's 
factor and compare to both the available experimental data and the predic- 
tions of the local density approximation, based on nuclear matter estimates. 
The agreement with the experiments is excellent, whereas the local density 
approximation is shown to sizably overestimate the Levinger's factor in the 

region of the medium nuclei. 

FACS number(s): 21.65+f; 21.60.Gx; 13.75.Cs 
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I. INTRODUCTION 



Within the Levinger's quasideuteron (QD) model [1-3] the nuclear photoabsorption cross 
section cr^Ej), above the giant dipole resonance and below pion threshold, is assumed to 
be proportional to the break-up cross section of a deuteron in hadronic matter, <tqd(E' 7 ): 



where E 1 is the photon energy and Vv is interpreted as the effective number of the nucleon- 
nucleon (NN) pairs of the QD type (see [4] and references therein ). Vv is written in the 
form 



where A and Z are the mass and atomic numbers of the nucleus and L is the so called 
Levinger's factor. Vv can be calculated for a given nuclear ground state wave function, thus 
allowing for a microscopic interpretation of the phenomenological Levinger's factor. 

The value of L has been extracted from experiments according to the following two 
models: i) the Levinger's model [5], in which o-qd(E^) is taken as the deuteron cross section 
damped by an exponential function, taking care of Pauli blocking of the final states available 
to the nucleon ejected from the QD, and ii) the Laget's model [6], which associates <tqd(£' 7 ) 
with the transition amplitudes of virtual (ir + p)-meson exchanges between the two nucleons 
of the QD pair. 

Both models provide satisfactory fits of photoreaction data in heavy nuclei, but yield 
different values of the Levinger's factor, Ll cv (^4) and LLa g ct(^4), LL agc t(^4) being ~ 20% 
larger than L Lcv (A). 

The effective number of deuteron-like pairs, as well as of three- and four-body structures, 
in spherical nuclei has been investigated within the shell model approach in Refs. [7,8]. In a 
recent paper [9] (referred to as I hereafter) we have analyzed the properties of deuteron-like 
structures in infinite symmetric nuclear matter (NM), described by a hamiltonian containing 
the realistic Urbana v u NN potential and the Urbana TNI many-body potential [10]. A 



o~a(E^) = V v o- QD (E 1 ) , 



(1) 



Vv = L 



Z(A - Z) 
A 



(2) 
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correlated wave function having spin-isospin dependent, central and tensor correlations has 
been used within the correlated basis functions (CBF) theory to compute the QD distribution 
function in matter and extract the NM Levinger's factor at equilibrium density, L^m = 
11.63, to be compared to the empirical estimate, L expt (A = oo) = 9.26. 

CBF theory has established itself as one of the most effective tools to realistically study, 
from a microscopic viewpoint, properties of infinite matter of nucleons ranging from the 
equation of state [11,12] to the momentum distribution [13] and the one- and two-body 
Green's functions [14-17]. In the last decade these studies have been successfully extended 
to deal with finite nuclei [18-22]. 

In this paper we extend the CBF many body approach, used in I for NM, to evaluate ab 
initio the momentum distribution, P D (k D ), and the total number per particle, Vt>/A, of QD 
pairs in the doubly closed shell nuclei 16 and 40 Ca, described in the Is coupling scheme. 
From Vt>/A we then extract the corresponding Levinger's factors. 

In Section II we review and generalize the CBF approach to the QD distribution in 
terms of the overlap between the nuclear and deuteron ground state wave functions. In 
Section III we compute the QD distribution and Vv in nuclei described by a realistic hamil- 
tonian including the modern Argonne v' 8 [23] and the Urbana IX [24] models of two- and 
three-nucleon potentials, respectively. The correlated nuclear wave function contains cen- 
tral and tensor correlations, as in Ref. [21]. The results are compared with the analogous 
NM quantities, obtained in I. We also evaluate the Levinger's factors, and compare them to 
the experimental values, as well as to those derived using the local density approximation 
(LDA) and the NM results of I. Summary and conclusions are given in Section. IV. 

II. QUASIDEUTERON DISTRIBUTION 

Following the approach developed in I, in a A-nucleon system the distribution of QD 
pairs whose center of mass is in the orbital state specified by the quantum number X can 
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be written 

Pd{X) = —^—{A\(a a D ) ] (X)a a D (X)\A) , (3) 

IdD + -L 

where \A) denotes the A-body ground state and J D = 1 is the spin of the deuteron. The 
operator a(af)%(X) annihililates (creates) a deuteron with the quantum number X in the 
a — 1, 2, 3 cartesian state. By introducing a complete set of intermediate (A — 2)- particle 
states and exploiting the completeness relation J2 n \n(A — 2)){n(A — 2)| = 1, we can recast 
Eq.(3) in the form 

P D (X) = J2(A\(a a D y(X)HA - 2))(n(A - 2)\a%(X)\A) . (4) 

In configuration space the above expression takes the form: 

J dRdh 1 d\ 2 d\ v dh 2 ^\{v 1 ,T 2 ,R)^ a D {X ] v l ,v 2 ) (y«(X;rr,r 2 ,))*y A (r v ,T 2/ ,R) , 

where R = (r 3 , . . . , r^), is the normalized nuclear ground state wave function and \&^) i s 
the deuteron wave function (DWF). 

The DWF can be split into its center of mass and relative motion parts according to 

V%(X; r u rj) = * Acro (I; Ry) VS,«i(v)|00> , (6) 

where = (r^ + r 3 -)/2, = — r^, 1 00) is the spin-isospin singlet NN state and 



u D {r l3 )o? - ^^T^>f 



(7) 



Ud{t) and Wd{t) being the t = and £ = 2 components of the deuteron wave function, 
whose normalization is given by 

POO 

/ r 2 dr [u 2 D (r) + w 2 D (r)] = 1 , (8) 
Jo 

In Eq.(7) cr° are the spin Pauli matrices, while the tensor operator reads 

r*(ry) = 3r*rg - ^ . (9) 
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Using the above definitions Pd{X) can finally be rewritten as: 

- ^j— (10) 

J dhidh^dh-rd 3 ^ D:Cm (X; R 12 )pg ) (n> r 2 ; r v , r 2 ,)V* Dcm (X; R V2 >) , 
where p$(ri, r 2 ; rv, r 2 /) is a generalized two-body density matrix defined by 

p%\r 1 ,T 2] r 1 ,,r 2 ,) = A(A-l) (11) 



dit!^(r 1 ,r 2 , J R)^(12)|00)(00|(^(l'2')) *A(r v ,r v ,R), 

where summation over the repeated indices is understood. 

The sum over X yields the total number of QD pairs in the nucleus, Vj>, thus allowing 
for a direct estimate of the Levinger's factor, L, to be compared to the empirical values 
resulting from phenomenological analyses [25,26] of photoreaction data [27,28]. 

A realistic A-body wave function, accounting for both short- and intermediate-range 
correlations induced by the strong nuclear interaction, is given in CBF theory by 

9 A (R)=s\Y[F(ijj\* (R), (12) 

i<j 

where R = (r 1; . . . , r^), S is a symmetrization operator and <£> is the Slater determinant of 
single particle orbitals (f> a (i), which are eigenfunctions of a suitable single particle hamilto- 
nian. For nuclear matter, the orbitals <f> a (i) are plane waves corresponding to a noninter- 
acting Fermi gas of nucleons with momenta |k| <kp — (G^Pnm/^) 1 ^ 3 , v = 4 and p NM are 
the NM spin-isospin degeneracy and density, respectively. 

The two-body correlation operator, F(ij), is given by the sum of 6 central and non- 
central spin-isospin dependent components, 

F(lj) = fMj) + faiUj)^ ■ <Tj) + f T {r lj ){T i ■ Tj) + far(r i] )((T l ■ (Tj)^ ■ T j) 

+f t (r ij )T a p(v ij )<j?o-P + ftrirjTap^ofafiTi ■ r 3 ) , (13) 

where the f p (r) correlation functions are variationally fixed by minimizing the ground state 
energy [21,29,30]. All the correlation functions heal to zero, except f c (r — > oo) — > 1. 



The generalized two-body density matrix p$ can be expanded in a series of terms in- 
volving an increasing number of nucleons by means of cluster expansion techniques [31]. In 
I the dressed leading order approximation (corresponding to the cluster diagram shown in 
Fig. 1 of I) was used to evaluate the momentum distribution of QD pairs in nuclear matter. 
The validity of this approximation has been satisfactorily checked in CBF calculations of 
the NM responses [32,33] and Green's functions [14,15]. In Ref. [22] the one-body density 
matrix of the N = Z doubly closed shell nuclei 16 O and 40 Ca has been computed using 
the correlation operator of Eq. (13), the realistic Argonne -Ug+Urbana IX interaction and 
the Fermi hypernetted chain/single operator chain (FHNC/SOC) diagrams resummation 
method [29,30]. Here we extend the approximation employed to calculate p^ in I to these 
two nuclei. 

f<y\ 

In the dressed leading order approximation p D (ri, r 2 ; iv, r 2 ') is given by 

pg^r^rv,^) » 2 J^l p m{ ruVv )i:{v 12 ,v V2 ,)p^{v 2 ,v 2l ), (14) 
where p^(r 1: ri>) is the one-body density matrix [22] and 



E(r 12 , v V2 ,) = -Tr Ft(l'2')< rd (l / 2 / )n o^^(12)F(12) (1 - P.P T ) ■ (15) 



n 00 and P a P T are the projector onto the (ST) = (00) two-nucleon state and the spin-isospin 
exchange operator, respectively (see I for details). 

By explicitly evaluating the trace in Eq. (15) in spin-isospin saturated systems, one gets 

S(r, r') = i; [U(r)U(r') + W(r)W(r')Q(r ■ ?)] , (16) 

with Q(x) = (3x 2 — l)/2, U(r) = u D (r) — Au(r) and W(r) = wd(t) — Aw(r). The Au(r) 
and Aw(r) functions account for the medium correlations effect on the bare components of 
the DWF. Their explicit expressions, in terms of the correlation functions, are given in I. 

Similarly to what is done to obtain the one-body momentum distribution in a nucleus, 
we consider the cm. orbital to be a plane wave with momentum ko in a periodical box of 
volume Q, 

^ Acm (k D ; Rij) = ^_ . (17) 
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As a consequence, for the QD momentum distribution (MD) we get: 



V v (k D ) = QP D (k D ) (18) 
= 7TT- / ^nd^rrd 3 ^ J*nin 12 -n ll2 ,) p (i) ^ ri ,) S ( ri2> ri , 2 ,)p (1) (r2, r 2 ,) . 

2 47T y 

This expression reduces to Eq. (13) of I in nuclear matter. Note that, in principle, different 
basis functions for the cm. orbitals, describing the spatial distribution of deuteron-like 
clusters inside the nucleus, can be used. 

In order to evaluate Vv(k D ) we define the function iV(ri,iv) through the relation [22] 

^(n, nO = n{t u nO £ xUi)x^(i') , (19) 

(JT 

where Xo-r(l) is the spin-isospin single particle wave function. Vx>(k D ) can be written in 
terms of the Fourier transforms of N, U and W as 

VvOcd) = y q ^ / dW (20) 

* (¥ ~ k> ¥ + k ') + w ( k ) w ( k ')Q& ■ ^ + k ^ - k ') • 

N(k, k') is related to JV(r, r') through 

iV(r, r') = / t^jfed 3 ^ e - l ( k r - k ' r ') N(k, k') , (21) 

and £7(/c) and J¥(/c) are given in I. 

In spherically symmetric nuclei spin and isospin indices are saturated and N(k, k') can 
be expressed in terms of Fourier-like transforms of the natural orbits (NO), <Pni°( k ) [22]: 

O/ i 1 ^ ^ 

N <& k ') = £ ^r p <( k • k ') n ^ ' ( 22 ) 

where denotes the /-th Legendre polynomial and <Pni°( k ) * s related to the configuration 
space NO, 0^°(r), through 

<ft°(k) = (2ir)-"* JdPrMkr)<^°(r), (23) 

ji(kr) being the spherical Bessel functions of order /. 
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The NO and their occupation numbers, n n [, are obtained by first expanding the one-body 
density matrix in multipoles, 



P (1) (ri,iv 



47T 



P ; (^-r)pr(n,nO, 



(24) 



and then diagonalizing p|^(ri, r±>), 



P?\r 1 ,r v ) = vY J n n l €?{r 1 )<p N n °{r v ) 



(25) 



The NO normalization is 



(26) 



In the independent particle model (IPM), ^a(R) = &o(R), and n^f M = 1, 4>nF = ^nz f° r 
occupied states, whereas n^f M = for unoccupied states. Deviations from IPM provide a 
measure of correlation effects, as they allow higher NO to become populated with n n i ^ 0. 
Using Eq. (22) we obtain: 



VvQld) = ^2 UaUa ' V v a \^o) , 

a,a' 



with 



16 4tt 



s=-2 



where a = (nlm), 



2 



+ k] /7(A;)0r 



k , 



and 



0X(q) = C u (g)>Wq), 



(27) 



(28) 



(29) 
(30) 

(31) 



^Zm(q) being the spherical harmonics. 



8 



III. RESULTS 



Last generation NN potentials are able to fit deuteron properties and the Nijmegen 93 
nucleon-nucleon scattering phase-shifts [36] up to the pion-production threshold (~ 4000 
data points) with a x 2 ~ 1- The Argonne v±$, belonging to this generation, is given by the 
sum of 14 isoscalar and 4 isovector terms, including charge- symmetry and charge-invariance 
breaking components [23]. In this work we have used a simpler NN potential, referred to as 
Argonne v' 8 , obtained from the the full Argonne v\$ retaining only the first eight operatorial 
terms, corresponding to those shown in Eq. (13) plus spin-orbit and spin-orbit/isospin. The 
Argonne v' 8 is constructed in such a way to reproduce the isoscalar part of the full vi$ in 
the S, P and 3 Di waves and the 3 D 1 - 3 Si coupling. The v' 8 parameterization, while allowing 
for a fully realistic NN interaction, makes the use of modern many-body methods, like CBF 
[21,34] or quantum Monte Carlo simulations [24,35] much more practical. It has been found 
that the differences between Argonne v' 8 and the full v 18 contribute very little to the binding 
energy of light nuclei and nuclear matter, and can be safely estimated either by perturbation 
theory or from FHNC/SOC calculations. 

It is well known that, to quantitatively describe the properties of nuclei with A > 2, 
modern NN interactions need to be supplemented with three-body forces. The Urbana IX 
(UIX) model provides a very good description of the energies of both the ground and the 
low-lying excited states of light nuclei (A < 8). In the present calculations we use Argonne 
v 8 + UIX interaction, which will be referred to as as the AU8' model. This interaction has 
already been used in the variational FHNC/SOC calculations of Ref. [22] as well as in the 
quantum Monte Carlo simulations of Ref. [24] . 

For the single particle wave functions, (j) a (i), entering the shell model wave function $ , 
we have solved the single particle Schrodinger equation with a Woods-Saxon potential, 

Vws(r) = — V ° • (32) 

1 + exp [(r - R )/a \ 

In principle, the parameters of the correlation functions, f p (r), and of the Woods-Saxon 



9 



potential may be both fixed by minimizing the ground state energy. This complete mini- 
mization was performed for the AU8' model in Ref. [21], and provided a binding energy per 
nucleon of B/A = 5.48 MeV in 16 and B / A = 6.97 MeV in 40 Ca (the experimental values 
are 7.97 MeV in 16 and 8.55 MeV in 40 Ca). These differences are compatible with the 
results of nuclear matter calculations at saturation density, p^M = 0.16 fm~ 3 , carried out 
with the same hamiltonian. In fact, the FHNC/SOC nuclear matter energy per nucleon is 
Enm/A = —10.9 MeV [21], to be compared to the empirical value of —16 MeV. 

However, the calculated root mean square radii of the two nuclei turned out to be R = 
2.83 fm in 16 and R = 3.66 fm in 40 Ca, showing a difference of ~ 5% with the experimental 
values, -Rexpt = 2.73 fm and -R ex pt = 3.48 fm, respectively. Moreover, the one-body densities 
were not in close agreement with the experimental ones. In order to take care of this feature 
of the variational approach, a set of single particle wave functions providing an accurate 
description of the empirical densities was chosen, and the energy was then minimized with 
respect to the correlation functions only. The resulting radii were R = 2.67 fm ( 16 0) and 
R = 3.39 fm ( 40 Ca), with a density description very much improved. The energies obtained 
by this partial minimization procedure were B/A = 5.41 MeV in 16 and B/A = 6.64 MeV 
in 40 Ca, largely within the accuracy of the FHNC/SOC scheme. Here, we adopt this same 
wave function, whose parameters are given in Table V of Ref. [21]. 

The structure of the NO in 16 and 40 Ca is discussed at length Ref. [22]. Here we limit 
ourselves to recall some of their main characteristics. The effect of correlations is mostly 
visible in the Is orbital, where the NO are larger than the shell model ones at short distances, 
resulting in stronger localization. The influence on the shape of the other occupied shell 
model orbitals is negligible. The occupation of the NO corresponding to the fully occupied 
shell model states is depleted by 9.6% in 16 and by 10.5% in 40 Ca, with a maximum 
depletion of ~ 22% for the 2s state in 40 Ca. As a consequence, the lowest mean field 
unoccupied states become sizably populated {ri2 S ~ n 2p ~ ri\d ~ 0.02 in 16 and n 3s 0.05, 
n 2p ~ 0.02, n 2 d ~ 0.03 in 40 Ca). These two effects are largely due to the presence of the 
tensor correlation. 
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Fig. 1 shows the behavior of U(r) and W{r) in 16 0, 40 Ca and nuclear matter, evaluated 
using the hamiltonian AU8'. For comparison, we also show the bare components of the 
Argonne v' 8 DWF. It appears that the main differences between deuteron and QD occur at 
r < 2 fm. At small relative distances both U(r) (r < 1 fm) and W(r) (r < 0.5 fm) are 
slightly suppressed with respect to ud{t) and wo(r). On the contrary, they are appreciably 
enhanced at larger distances. These effects are more visible for the lightest nucleus. 

The differences between nuclear matter and nuclei mostly disappear in the Fourier trans- 
forms, \U(k)\, \W(k)\, \uo(k)\ and \wo(k)\, whose behavior is displayed in Fig. 2. The nu- 
clear medium shifts the second minimum of \uo{k)\ towards lower values of k, as obtained 
in I for nuclear matter with the Urbana V\± potential. The Argonne v' 8 \wr>(k)\ does not 
exihibit any diffraction minimum, which, however, appears in |W(fc)|. 

The distribution of deuteron pairs with total momentum k D , VvQc-d), resulting from our 
approach is displayed by the solid line in Fig. 3 for 16 O and in Fig. 4 for 40 Ca. The following 
comments are in order: 

(i) NN correlations introduce high momentum components in the distribution. The full 
V-dO^d) is strongly enhanced with respect to V^ PM (k D ) at large |ko|, and it is cor- 
respondingly depleted at small |ko|. The depletion is mostly due to the non-central 
tensor correlations. 

(ii) The effect of state-dependent correlations is large, as one can see by comparing the 
full Vvik-D) with the Jastrow model V^ik^) (obtained by retaining only the scalar 
component in the two-body correlation operator (13)). 

(iii) The tail of V-dO&d) is appreciably different from that of nuclear matter. At |kc>| = 4/cp 
the difference is still a factor ~ 10 for both 16 and 40 Ca. 

Fig. 5 displays the convergence of Vviko) in the number of natural orbits included in 
the sum of Eq. (27). Full convergence is reached with the inclusion of orbits up to 5/ for 
both 16 and 40 Ca. The figure shows that, in the case of 40 Ca, the tail of Vi>(k D ) is still 
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~ 10 times too small if only orbitals up to 3d are included. 

The contributions to Vx>(k D ) coming from the various orbitals are displayed in Fig. 6 for 
16 and on Fig. 7 for 40 Ca. They are also compared with the corresponding results obtained 
within the IPM and Jastrow models. The effects of state-dependent correlations is large in 
all channels, and particularly in the highest ones. 

The total number of pairs of the QD type in both finite nuclei and nuclear matter, Vt> 
is obtained by integration of Vv^d) over k D : 



where the factor 3 in the r.h.s corresponds to the spin multiplicity of the deuteron, 2 Jd + 1. 

We have repeated the calculations for nuclear matter by using the AU8' interaction. 
The result Vv(NM)/A = 2.707 (the corresponding Fermi gas model result is 3.382) should 
be compared to the value 2.895, obtained in I with the Urbana vu two-nucleon plus the 
Urbana TNI many-body forces [10] (which will referred to as the UU14 model). The corre- 
sponding numbers for 16 and 40 Ca turn out to be much smaller: Vt>{ w O)/A = 1.090 and 
V v rCa)/A = 1.370, respectively. 

The Levinger factor is easily obtained from Vv/A by means of Eq. (2). As we are dealing 
with symmetric matter (N = Z = A/2), L(A) = AVv/A. Our estimates L f6 ( 16 0), L f6 ( 40 Ca) 
and Lf 6 (NM) for 16 0, 40 Ca and nuclear matter, corresponding to the / 6 correlation model, 
are reported in Figs. 3 and 4. These results are not too different from the values obtained 
within the independent particle and Jastrow models. This fact actually implies that the high 
momentum tail of Vt>(^d) is not relevant for the calculation of the Levinger factor L. It has 
to be stressed that the Jastrow model turns out to consistently underestimate the Levinger 
factor. The spatial structure of np pairs having the deuteron quantum numbers has been 
investigated in Ref. [37] in light (A=3,4,6 and 7) nuclei and 16 using a variational Monte 
Carlo approach and the Argonne v is two-nucleon and Urbana IX three-nucleon potentials. 
The estimated Levinger factor for 16 is L V mc( 16 0) = 4.70, comfortably close to our value, 




/ 



d 3 k D Vy(k D ) 
(2vr) 3 A 



(33) 



L f6 ( 16 0) = 4.36. 
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Our results for the Levinger factors are summarized in Fig. 8, where they are also com- 
pared with the available experimental estimates. The agreement with the photoreaction data 
of Ahrens et al. [28] for the case of 16 and 40 Ca is rather impressive. The "experimental" 
value, Lexpt(oo) = 9.26, deduced from the phenomenological formula 

L cxpt (A) = 13.82— ^-3- , (34) 
R 6 [tin J 

reported in Ref. [25], is ~ 15% smaller than our theoretical value. In I, the surface con- 
tribution to L(A) has been estimated exploiting the calculated enhancement factor of the 
electric dipole sum rule for finite nuclei, )C, [38], obtained using CBF theory and LDA. The 
enhancement factor is related to experimental data on photoreactions through the equation: 

1 + /C cxpt = — / a A {E 1 )dE 1 , (35) 
°o Jo 

where a = 60 [Z(A — Z)/A] MeV mb and is the 7r-meson production threshold. By 
using the same parameterization as in I for the surface term, we get: 

L LDA (A) = 10.83 - 9.76 A~ 1/3 , (36) 

for the AU8' interaction. L LD \(A) is displayed on Fig. 8. LDA turns out to be not satis- 
factory for medium nuclei, such as 16 and 40 Ca. Fig. 8 also report L Lcv (^4) and L Lagct (A), 
as extracted [25,26] from the available experimental data on photoreactions. The computed 
Levinger's factors are almost A-independent for heavy nuclei (A > 100), and result to be 
~ 15% larger than L Lev (A) and ~ 25% smaller than L La gct(^4)- Such disagreement between 
theory and experiment is likely to be ascribed to the sizable tail contributions to the electric 
dipole sum rule, absent in the definition of Eq. (35). 

IV. CONCLUSION 

The Correlated Basis Function theory of the two-body density matrix has been applied to 
microscopically compute the distribution of QD pairs carrying total momentum ko, VvQc-d), 
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in doubly closed shell nuclei 16 and 40 Ca and nuclear matter, starting from the realistic 
Argonne v' 8 plus Urbana IX potential. 

It has been found that NN correlations produce a high momentum tail in VvO^d) and, 
correspondingly a depletion at small kp for both nuclei and nuclear matter. These effects 
are mainly due to the presence of the state-dependent correlations associated with the 
tensor component of the one pion exchange interaction. Contrary to what happens for the 
one-nucleon momentum distibution, the tail of VvQc-d) sizably differs from that of nuclear 
matter. 

Summation of Vv^d) over ko provides the total number Vv of QD pairs, and, conse- 
quently, allows for an ah initio calculation of the Levinger's factor ~L(A). The CBF result 
for nuclear matter is significantly reduced with respect to the value obtained in I with the 
Urbana V\± plus the Urbana TNI many-body forces. The corresponding Levinger factors for 
16 and 40 Ca are much smaller than the nuclear matter value and in very good agreement 
with the available photoreaction data analyzed within the quasideuteron phenomenology. In 
addition, our results show that LDA overestimates L(A) in the region of the light-medium 
nuclei. 

The ~L(A) resulting from the full calculation are relatively close to the corresponding 
values obtained within the IPM and Jastrow models. Actually, the high momentum tail of 
V-dO*-d) gives a small contributions to the Levinger factor. This feature indicates that the 
approximation used in our calculation (which amounts to including only diagrams at the 
dressed lowest order of the FHNC cluster expansion) is fully adequate. However, it should 
be noticed that the Jastrow model underestimates the Levinger factor. 

In addition, the analysis described in this paper shows that when a deuteron is embedded 
in a nucleus, or in nuclear matter at equilibrium density, its wave function gets appreciably 
modified by the surrounding medium. While in the case of the S'-wave component the 
difference is mostly visible at small relative distance (r < 1 fm), the D-wave component of 
the QD appears to be significantly enhanced, with respect to the deuteron W£>(r), over the 
range r < 2 fm. This effect is particularly evident in the lightest nucleus. 
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FIGURES 




r [fm] 

FIG. 1. Radial components U(r) and W(r) of the AU8' QD wave functions in 16 0, 40 Ca and 
nuclear matter. Upper panel: the solid and dashed lines show the radial dependence of Ua(t) for 
16 and 40 Ca, respectively. The dot-dashed and dotted lines correspond to the nuclear matter 
UNM( r ) and the bare ii£)(r). Lower panel: as in the upper panel for the d-wave components of 
the QD and deuteron wave functions. 
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FIG. 2. As in Fig. 1 in momentum space. 
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1^1/(2**) 

FIG. 3. Momentum distribution of QD pairs in 16 O as a function of the total momentum 
|ko| (see Eq. (27)). The solid, dashed and dash-dotted lines are the results obtained within the 
/g and Jastrow correlation models and IPM, respectively. The short-dashed line displays the 
momentum distribution of the QD in nuclear matter at equilibrium density, pnm = 0.16 fm~ 3 . 
The insert shows a blow up of the region |k£>|/(2fcp) < 1, plotted in linear scale. The Levinger 
factors, ~L(A), for the various calculations are also reported. 
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FIG. 4. As in Fig. 3 for 40 Ca. 
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FIG. 5. Convergence of Vt>0<-d)/A for 16 and 40 Ca in the number of natural orbits included 
in the summation of Eq. (27). The results have been obtained within the f$ correlation model. 
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FIG. 6. Diagonal contributions from the [Is — 2d] orbitals to Vt>0<-d)/A for 16 0. The solid, 
dashed and dash-dotted lines refer to the Jq, Jastrow and IP models, respectively. The occupation 
numbers n„i( 16 0) computed with the f$ correlation [22] are also reported. The numbers in square 
brackets refer to the Jastrow case. 
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FIG. 7. As in Fig. 6 for 40 Ca. 
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FIG. 8. Levinger's factor L(A) for 16 0, 40 Ca and nuclear matter (shown by the arrows for the 
UU14 and AU8' forces). The filled circles, the empty circles and the triangles show the Levinger's 
factors obtained within the f§ and Jastrow correlation models and the IPM, respectively. The LDA, 
as discussed in the text, is also reported (solid line). The phenomenological values of LL ev (^) 
corresponding to the photoreaction data of Lepretre et al. [27] (squares) and Ahrens et al. [28] 
(crosses and diamonds) are taken from Ref. [25]. The empirical values of LL age t (A) , represented by 
circles in the heavy nuclei region, are from Ref. [26]. 
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